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ABSTRACT 

Motivated by the detection of a recent outburst of the massive luminous blue variable 
LMC-R71, which reached an absolute magnitude My = —9.3 mag, we undertook a system- 
atic study of the o ptical variability of 12 68 massive stars in the Large Magellanic C loud, using 
a rec ent catalog by iBonanos et al. I (l2009h as the input. The ASAS All Star Catalog |Po imahskil 



2002) provided well-sampled light curves of these bright stars spanning 10 years. Combining the 



two catalogs resulted in 599 matches, on which we performed a variability search. We identified 
117 variable stars, 38 of which were not known before, despite their brightness and large ampli- 
tude of variation. We found 13 periodic stars that we classify as eclipsing binary (EB) stars, eight 
of which are newly discovered bright, massive eclipsing binaries composed of OB type stars. The 
remaining 104 variables are either semi- or non-periodic, the majority (85) being red supergiants. 
Most (26) of the newly discovered variables in this category are also red supergiants with only 
three B and four O stars. 



Subject headings: Magellanic Clouds — binaries: eclipsing — stars: massive — stars: variables: general 
1. INTRODUCTION 



1 Hubble and Carnegie-Princeton Fellow 

Studies of stellar variability in the Magellanic 
Clouds have a long and venerable history, start- 
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ing with a 1904 Harvard College Observatory Cir- 
cular by E. C. Pickering stating that "The two 
Magellanic Clouds have long been objects of care- 
ful study, on account of the extraordinary physi- 
cal conditions which prevail in them. They have 
not, however, heretofore been known as regions 
in wh ich variable stars are numerous" ((Pickering 
19041 ). That circular, although formally authored 
by Pickering, in fact reported 152 new variable 
stars in the Large Magellanic Cloud found by 
Henrietta Leavitt, work that was soon followed 
by "843 New Variable Stars in the Sma ll Mag- 
ellanic Cloud" ([Leavitt fc Pickering! Il905l). "1777 
Variables in the Magellanic Clouds" (|Leavitdl 19081 ) 
and culminated in the formulation of the Cepheid 
period-luminosity relatio n (the "Leavitt Law", 
Leavitt fc Pickering|l912 ) , a result of fundamental 
and continued importance in astron o my and cos- 



Clouds and supe rnovae light echoes in the LMC 
(jRest et al.ll2005! ) . Some microlensing events were 
also found in the Clouds, but clearly not enough 
to populate the halo of our Galaxy with MA- 
CHO s (e.g.. lAlcock et all 1200 " 
2007t IWvrzvkowski et al.ll200 " 



Tisserand et al 



2001; Macrietal 



molo gy (e.g., iFreedman et al 
120061) . 

Later studies of stellar variability in the Clouds 



Given all that tremendous effort and many hun- 
dreds of papers written about stellar variability 
in the Magellanic Clouds, is there still room for 
more discoveries? As is almost always the case 
in astronomy, the answer is "yes". One path 
is simply to go deeper, and some impressive re- 
sults have been reported recently, such as the de- 
tection of_extragalactic_W UMa binaries in the 
LMC (jKaluzny et al.l 120061 ) using ground-based 
data from a large aperture telescope. Another 
path is to find variable stars in the densest parts 
of the Clouds, which are now accessible using 



the H ubble Space Telescope (e.g.. iFiorentino et al 



2008). Still yet another path is to take advantage 



of the fact that we have been observing stars in 



included many papers on Cepheid s fe.g.. lCaldwell fc 
19861: iMadore fc Freedmanlll99l|) . but there were 
also studies of long peri od variables (e.g. JWood et al 
19831 [Feast et al.lll989l) . irregular variables of the 
S D oradus type (lumin ous blue variables, LBVs; 
e.g.. IStahl et al. 19831) and R R Lyr stars (e.g., 
Thackeravf 119581 : IWalkeJll992l) . just to name the 



most popular subjects. By virtue of being so close 
to us and being relatively unobscured, the Magel- 
lanic Clouds continue to provide an ideal terrain 
for stellar variability studies. 

The idea of IPaczvhski ( 19861 ) to use the Mag- 
ellanic Clouds as targets for microlensing surveys 
to detect dark matter in the form of MACHOs 
has led to a new and glorious era in variability 
studies in general, not just in the Clouds, and 
has resulted in a wide variety and large number 
of papers reporting on variability in the Large 
Magellanic Cloud (LMC) and the Small Magel- 
lanic Cloud (SMC), includ i ng work on Cepheids 



Coulschble Magellanic Cloud s for a long time now (e.g., 
iKozlowski et al.l[2 010 N ). even as long as a hundred 
years (jGrindlav et al.1 120091 ) . But the path taken 
in this paper, namely to go shallower, and to cover 
all of the Clouds (the SMC will be discussed in a 
future paper), again owes its existence to Bohdan 
Paczyhski, whose vision also led to the All Sky 
Automated Survey (ASAS) project led by Grze- 
gorz Pojmahski, who created it and has kept it 
functional for the past 10 years. 

A particular motivation for our study was pro- 
vided by the recently reported outburst of an LMC 
LBV R710, which during the last 5 years has 
brightened by more than 2 magnitudes to a vi- 
sual magnitude of abou t 9.0, becoming the b right- 



est star in the LMC dGamen et al.l l2009h . As 

it turns out, the ASA S (|Poimahski Il997t 12002 : 
Pojmahski et al. 20051 ) photometric database con- 



Alcocket all Il995t lUdalski et all 1 19991 



Soszvhsk i et al.l I2008D. eclipsing binaries (e.g., 
Alcock et al.lll997t IWvrzvkowski et al.ll2004l). RR 
Lyr s t ars (e.g-.lAlcock et al.lll996HSoszvhski et al 



varia bility (e.g. 
20011). 



20031 ISoszvhski et al.l 120091) all k i nds of stellar 



Dobrzvcki et al.l [2003 



Cook et al. 1995; Zcbruh et al. 


quasars (e.g., Geha et al. 


2003: 


2003: Koz 


owski et al.l 


2010). 


iGarg et al. 


20071) behind the 



tained a detailed F-band light curve of that star 
(see Figure [T]), obtained and reported in an au- 
tomatic fashion by ASAS. Motivated by the fact 
that the brightening of R71 started several years 
ago, but was first reported only in 2008, we de- 
cided to investigate light curves of other bright 
stars in the LMC, using the ASAS data. R127 
is another example of a ninth magnitude LBV in 
the LMC that underwent dramatic change start- 



1 R71 = ASAS 050207-7120.2 
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ing in 2005 that was recorded by ASAS^ (see Fig. 
H), and was not noti ced until two years later by 
Walborn et al.l ( 2008 ). who began monitoring it 
spectroscopically. There are bound to be other 
such examples of bright LMC variables already 
recorded in the AS AS database. We therefo re used 
the recent catalog of Bonanos et al.l ( 20091 ) as an 
input catalog to characterize the variability prop- 
erties of these known LMC massive stars. Con- 
firming our suspicion, there are not only many 
spectacular light curves among known variables, 
but also a significant fraction of new interesting 
objects. 

This paper has the following structure. In Sec- 
tion 2, we briefly describe both the input catalogs 
and the matching procedure. Section 3 introduces 
our variable star selection methods and discusses 
the identified variables. In Section 4, we present 
color-magnitude diagrams (CMDs) of the variable 
stars and discuss special cases. We summarize this 
work in Section 5. 

2. DATA 

2.1. Input Catalogs 



We use the catalog of iBonanos et al. (l2009h . 



hereafter B09, as the source of bright, massive ob- 
jects associated with the LMC. It contains pho- 
tometric measurements for 1268 massive stars in 
several visual and infrared wavelengths ( UB VI, 
JHKs, and 3.6-24 /an), as well as their spectral 
classifications. The area covered by this catalog is 
delimited by -72° < S < -64° and 70° < a < 

90°. 

The ASAS^ (|Pojmahskill2002l ) has been moni- 
toring the southern sky (S < +28°) since 2000 and 
by now the average number of data points per star 
in the LMC region is 550 in the V band and 350 
in the I band. The V magnitude range covered 
by the ASAS telescopes is around 8 - 14.5, which 
ideally suits a study of the brightest objects in the 
LMC. The number of objects visible in the LMC 
direction, within the coordinates listed above, is 
about 28 000, many of them being Galactic field 
stars. 



2 R127 = ASAS 053644-6929.8 (although blended) 
3 http:/ /www. astrouw.edu.pl/asas 



2.2. Matching Procedure 

The main difficulty in combining the two cata- 
logs resulted from the large size of the ASAS pixel, 
which is ~14'/9. Thus, even though coordinates 
of obj ects in the ASAS catalogs are accurate to 
- 3" (|Poimanskilll997h and the positions within 
the B09 catalog are typically accurate to less than 
r~> 1", in very crowded regions one ASAS object 
may be composed of more than one (bright) LMC 
star. 

For this reason we performed catalog matching 
in two steps. First, each of the 1268 stars from the 
B09 catalog was assigned one closest counterpart 
from the ASAS catalog (28,000 sources), without 
putting any constraints on counterpart separation. 
Among 1268 matching entries, 761 turned out to 
be unique, which means that an ASAS star in such 
a pair was assigned only one counterpart from the 
B09 catalog. Within the remaining 507 match- 
ing entries, one ASAS star was assigned to either 
two or three B09 stars and as a result these 507 
matching entries contained only 162 unique ASAS 
sources. This was the case in the densest regions 
of the LMC, where due to low ASAS resolution ob- 
jects from the B09 catalog could not be resolved. 

In the second step we took the 923 (761 + 
162) ASAS stars selected in the previous step and 
checked how many B09 neighbors each ASAS star 
had within a radius of 15" (1), between radii 15" 
and 30" (2), and between 30" and 45" (3). Given 
that the aperture radius adopted by us is in most 
cases (~70%) 1 ASAS pixel, we assumed that 
the match could be considered unambiguous when 
there was only one neighbor within (1) and none 
within (2). There were 675 objects fulfilling this 
criterion. However, for about 30% of ASAS stars 
we have adopted apertures greater than 1 pixel, 
which could result in neighboring stars overlap- 
ping each other, so we decided to apply a more 
conservative criterion, that no stars were allowed 
in either zone (2) or (3). The number of matches 
was reduced to 599, which we accepted as our final 
sample. The choice of the aperture size described 
above is automatic and depends only on the mag- 
nitude of an object — for the faintest stars (13 - 
14 mag) we use the smallest aperture radius of 1 
pixel, while for the brightest (8 mag) we use the 
largest aperture radius of 3 pixels. 

It is important to stress here that the match has 
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been made only against the B09 catalog, which 
means that the foreground field stars are not ac- 
counted for. At the same time, the B09 catalog 
itself is far from a complete list of massive stars in 
the LMC, as it contains only the stars for which 
spectra exist. Therefore, we should expect some 
crowding problems even in our final sample. This 
issue is addressed in Section 3.3. 

The spatial distribution of the 599 matched ob- 
jects in equatorial coordinates is presented in Fig- 
ure [2] with black filled and open circles. Gray dots 
represent all ASAS objects in this region brighter 
than V = 14 mag (-19,000 of -28,000 points). 
The sample is dominated by foreground stars. 

3. RESULTS OF VARIABLE STAR SEARCH 

As already summarized in the Introduction, the 
LMC has been the target of variability searches 
since the beginning of the last century, with an 
increasing intensity over the last two decades. Re- 
cent catalogs of variable stars were mainly gener- 
ated using the data acquired by microlensing sur- 
veys such as OGLE and MACHO. However, these 
catalogs are complete for stars fainter than 13-14 
V magnitude, and are missing the brightest ob- 
jects. This oversight can be corrected with the aid 
of ASAS data. 

3.1. Variable Star Selection 

In order to investigate variability of the 599 
matched stars, we extracted y-band light curves 
from the online ASAS database, keeping only 
points with quality flag "A". Then we cleaned 
the light curves by removing points that lie more 
than 3cr from a local linear model fitted to each 
light curve, where a is the standard deviation. 
All cleaned light curves were subject to variabil- 
ity analysi s using the analysis of varian ces (AoV) 
algorithm ( Schwarzenberg-Czernvll989l ) in the fre- 
quency range 0.0001 — 8 c d _1 . Each light curve 
was then folded with a period equivalent to the 
frequency with the highest peak in the power spec- 
trum and visually inspected to eliminate possible 
spurious detections such as 1 day aliases. During 
the visual inspection other frequencies with high 
signal were examined and the best period was cho- 
sen. 

Figure [3] presents the standard deviation cry 
versus y-band mean magnitude with all 599 



matched sources marked with gray dots and vari- 
able stars with larger circles. There are a number 
of stars with high cry values that were not clas- 
sified as variable. These objects have an increase 
in brightness around HJD 2452200 and 2452600 
which was due to strong defocusing of the ASAS 
telescopes. 

During the inspection process we identified 
117 variable stars, of which 55 were not listed 
as variable in the SimbacQ database. We then 
matched these 55 s tars against the GCVS catalog 
(|Samus et a l. 2009) using VizieR0 and found that 
another 24 were known to be variable, which alto- 
gether gives 79 known variable stars and 38 new 
discoveries. 

3.2. Periodic Variables 

Among our 117 detections there are 13 periodic 
variables; five of them arc known binary systems 
(although BAT99-32 was only known as a spec- 
troscopic binary and its eclipsing nature has not 
been reported before), which are as follows: 



045515-6711.4 
Sk 



BAT99-6 = Bre y 5 



-67 18 (IPrevot-Burnichon et al.1 119811 ): 
ISeggewiss et al. ( 199lh have identified this 
object as a binary system with a tentative 
period of 1.99 d ays, based on a very poor 



light curve. Later [Niemela et~aL ( 2001 ) have 
shown it to be "a multiple system consisting 
of at least two pairs of short period bina- 
ries." We confirm the 2 day period of the 
eclipsing pair. 

052223-7136.0 = BAT9 9-32: a long known 
Wolf- Rayet (W-R) star (|Westerlund fc Smith 



19641 ). which has been reporte d to be a sin- 



gle lined spectroscopic binary (lMoffatlll989h 
and an EB with a p eriod of —1.91 days 
(<Seggewiss et al.lll99ll) . The discovery from 



ASAS observations that it is also an eclips- 
ing system with a period P —3.8 days that is 
twice as long as that found in the literature, 
explains the low radial velocity amplitude 
and unusually short period previously re- 
ported. High resolution spectra of this sys- 
tem would be very valuable as they would 



4 http:/ /simbad.u-strasbg.fr 
5 http:/ /vizier. u-strasbg.fr 
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allow for a determination of its parameters, 
given that very few W-R stars have good 
measurements of masses and radii. 



05260 6-6710.9 = Sk -67 105: iNiemela fc Morrell 
(1986) found it to be a massive double lined 
spectroscopic binary (the firs t extragalactic 
binary to be studied), while lHaefner et al 



(|1994h discov ered it is a near con t act ec lips- 
ing system. lOstrov fc Lapasset ( 2003 ) de- 
rived component masses of 48.3 M and 
31.4M Q . 



• 052727-6712.0 = Sk -67 117 = HV 2543: a 
known semi-detached double-lined EB with 
masses 25.63 ± 0.7 Mp, and 15.63 ± 1.0 M w 



(|Niemela fc Bassinoll994llOstrov etllB oOO) 



053442-6931.6 = MACHQ81.8763.8: discov- 
ered as an EB bv lAlcock etaD (|l997h . it was 
shown to be a massive binary with compo- 
nent masses of 41 M Q and 27 M Q by lOstrov 



(|2001|) 



The y-band light curves of these binaries, span- 
ning 10 years of AS AS observations are plotted in 
Figure [4] and their parameters are presented in Ta- 
ble 1. In combination with previously published 
photometric data, the ASAS observations can be 
used to improve the precision of the binary orbital 
periods, by extending the baseline of the observa- 
tions. 

Parameters for the remaining 8 objects (blue 
pentagons in Figure [2]) are listed in Table 2, and 
their light curves are plotted in Figure [5] They are 
most probably newly discovered EB stars and we 
have already started follow up spectroscopic obser- 
vations in order to verify their binarity. All have 
OB spectral types, which is consistent with their 
short orbital period. Three of them had been rec- 
ognized as variable stars long ago, one was misclas- 
sified as a long period variabl e (052506-712 7.9 — 
Sk -71 29 = LMC V252&9) bv lHughesl (Il989h . and 
the other two had no type specification: 05 1230- 
6727.3 = BI 98 (iButler fc Wavmanl Il974l) a nd 
053435-6945.7 = Sk -69 194 (IWesterlundll96lh . 

The latter, namely 053435-6945.7 (Sk -69 194) 
has the longest period of all EB candidates (P = 



6 "Poorly studied irregular variables of early (O-A) spectral 
type." GCVS 



12.2 days). According to Massev et al. ( 2000l ) it is 
a B0 la + WN star. The filling-in of the minima 
in its light curve (Figure [5]) most probably results 
from the presence of a bright (V « 11 mag) fore- 
ground K type star HD 269770, distant by 33" in 
the SW direction. 

The remaining five EB candidates were not pre- 
viously identified as variable stars. Their ASAS 
identifications are as follows: 



0456 17-6959.6 = Sk -70 18 (B0 m lConti et al" 
1986b 



0449 51-6912.0 = Sk -69 9 (06.5 IIF lJaxon et al 



2001) 



• 052055-6527.3 = Sk -65 47 (04 If), whose 
spectrum is available in the l iterature (Fig- 



ure 7 of iMassev et~ai1 120051) . The He II 



A4686 and Ha lines are not well fit by the 
model, which may be due to wind inter- 
action in the binary. The object was re- 
cently reclassified to 04 I(n)f+p type by 



Walborn et ail (|2010l ). Since 04 If stars are 



quite rare, measuring the parameters for this 
system would be valuable. 



• 05260 1-6703.1 = Sk -67 102 (B2 IIL lConti et al 
1986r ) 



• 054711-6759.0 = Sk -67 270 (B0. 5 V; 



Conti et al.l ll986) 



3.3. Non-periodic Variables 

We identified 104 variable stars that are either 
semi-periodic or non-periodic. Among them are 
25 variables that were already listed in the ASAS 
Catalog of Variable Stars (ACVS) as MISC-type 
variable stars, but 9 of these are not listed as vari- 
able in the Simbad database. Altogether there are 
33 objects that are not marked as variable in Sim- 
bad, GCVS, or ASAS. All of them are bright and 
many have large light curve amplitudes, so it is 
noteworthy that their variability had not been dis- 
covered previously. The reason that some of them 
do not figure in ACVS is that the part of the cata- 
log containing the LMC region was made when the 
light curves cons isted of 40-50 observations only 
(|Pojmahskill2002h . during which there was no sig- 
nificant change in the light curves. 
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The light curves of some interesting objects that 
are for the first time reported as variable are plot- 
ted in Figure |6l while Figure [7] shows long term 
light curves of a few known interesting variables; 
the last two are solely ASAS discoveries. In Figure 
[8] we present light curve stamps of all 104 semi- or 
non-periodic stars. 

The parameters of all 117 variable stars, both 
periodic and non-periodic are listed in Table 
3, supplemented with extinction values from 
Peicha fc Stanek ( 20091 ) where available (101 ob- 



jects). The next to last column contains informa- 
tion on whether the star is a known variable (by 
providing its ID) or if it is a new discovery (empty 
field). 

The last column provides information about 
crowding problems mentioned in the fourth para- 
graph of Section 2.2. It contains the number of 
bright companions and their approximate angular 
distance (in arcsec) from the variable star. These 
objects were selected by examining the Digitized 
Sky Survey (DSS) images around each variable 
star, using the Goddard Sky Viei^\ utility. The 
DSS exposures reach 2-3 mag deeper than ASAS 
observations, making them ideal for resolving ob- 
jects otherwise indistinguishable by ASAS. 

The search showed that 47 out of the 117 vari- 
ables have fairly bright close companions. Usually 
the presence of such neighbors will be reflected 
in the light curves by a significant scatter in the 
observations, as in the case of an EB Sk -69 194 
already mentioned in Section 3.2. Another two ex- 
amples are RSG stars 053627-6923.9 (one bright 
star ~ 10" away) and 052749-6913.3 (1 very bright 
star ~ 20" away) presented in Figure 

Table 3 is available as a whole in the electronic 
version of this article. 

Using the classification of massive stars adopted 
in B09 we have divided the stars in our sample 
into the following groups, sorted by the number 
of objects in the group (the abbreviation is given 
in brackets). For each class we list the number 
of variables, including periodic variables, and the 
total number of stars in our matched sample: 

- early B type stars (earlyB): 10 / 238, 

- O type stars (O): 10 / 156, 

- red supergiants (RSGs): 85 / 97, 



- late B type stars (lateB): 2 / 51, 

- Wolf Rayet stars (WR): 3 / 35, 

- A, F, G supergiants (AFG): 1 / 10, 

- supergiant B[e] stars (sgBe): 3/8, 

- luminous blue variables (LBVs): 3/3, 

- Be/X-ray binaries (bexrb): /I. 

As expected, RSGs constitute the majority of all 
variable stars (85/117=73%), and they have the 
largest variability percentage among all classes 
other than the LBVs (by definition). It is also 
difficult to confirm that the 12 RSGs that we do 
not report as variable are truly not variable. Their 
light curves show quite high scatter, which could 
indicate lower amplitude variability than is de- 
tectable by ASAS. The percentage of variables in 
other classes varies from about 4% for B type stars 
to about 37% for sgBe stars, but the number of ob- 
jects in several groups is too small to perform any 
global statistics. 

We will discuss some of these variable stars in 
the following section. 

4. VARIABLE STAR PROPERTIES 
4.1. Color Magnitude Diagrams 

In the left panel of Figure M we plot a CMD, 
showing V versus V — I magnitudes of all final 
matches between ASAS and B09 catalogs. V and I 
are the mean magnitude values based on the ASAS 
light curves. Gray open squares represent all 599 
matches, while black dots stand for variable stars. 
EBs are marked with big open circles. The ma- 
jority of variables have red colors, as noted in the 
previous section. At the same time almost all red 
stars turn out to be variable, which su pports pre 



vious studies of RSG variability (e.g., iKiss et al 



7 http: / /skyview. gsfc.nasa.gov/cgi-bin/query. pi 



2006, and references therein). All eclipsing vari- 
ables are among the blue stars. The right panel 
of Figure |9] presents V amplitude vs V — I color 
for the same sample with the same designations. 
The amplitude range is large, from ~0.1 to ~1.8 
mag. Non-variable stars (gray squares) are plot- 
ted to show the scale of observational scatter in 
the dense regions of the LMC. For these, the am- 
plitude Ay corresponds to the scatter in the data 
points, being roughly the 3xov value. The ma- 
jority of objects with high A\ values are among 
the faint stars (see Figure [3|). 

The star with the largest amplitude value, 



6 



namely ASAS 051450-6727.4 = [SP77] 37-35 = HV 
916, is a known variable supergiant of type M3 la 
(its light curve is plotted in Figure |7|), while R71 
with its dramatic rise is only fifth in amplitude and 
is marked with a cross in the right panel of Figure 
H The lowest amplitude variable is ASAS 045647- 
6950.4 = R66 which has been recognized as S Dor 
/ a Cy g type variable bv lvan Genderen fc Sterkenl 
([2002) and will be discussed in Section 4.2. 

Figure [TU] is a CMD diagram of Spitzer/IRAC 
[3.6] versus [3.6]-[4.5] infrared magnitudes of the 
117 variables, taken from B09 catalog. This figure 
is analogous to the Figure 2 of B09 where different 
spectral types are marked with different symbols, 
but here we additionally mark eclipsing variables 
with large open circles. Again we see that the ma- 
jority of variable stars are among RSGs, occupying 
the region of [3.6] < 9 and [3.6]-[4.5] < 0.4 mag. 
All EBs are either O- or early B-type objects, ex- 
cept for one W-R star. 

4.2. Special Cases 

Three of the sgB[e] stars, marked with blue di- 
amonds in Figure | 10| were previously reported as 
variable by Ivan Genderen fc Sterken] (|l999l l2002h 
and exhibit low amplitude variability: S22 = HD 
34664 = ASAS 051353-6726.9 (~ 0.21 mag), R66 
= ASAS 045647-6950.4 (~ 0.08 mag) and R126 
= 053626-6922.9 (~ 0.16 mag). We find evidence 
for a long-term S Dor-like oscillation, on the or- 
der of_sjx_^ea£S J _jn_S22 J _jhn that reported 
by Ivan Genderen k, Sterken ( 19991 ). who conclude 
that this star "is not only a B[e] star, but also a 
very weak-active LBV" . 



In the case of R66 van Genderen fc Sterken 



(|2002[ ) observed "a long-term wave with a possi- 
ble timescale of hundreds of days and an ampli- 
tude of ^0^06" based on their observations from 
1989 to 1991. This is the lowest amplitude variable 
in our sample, with Ay = 0.08 being consistent 
with previous estimates. We observe a pulsation- 
type oscillation with a period of 224 days, which 
is roughly four t imes longer than the 55 day pe - 
riod observed bv Ivan Genderen k Sterken! (|2002l ). 
attributed to a Cyg-type oscillations. We do not 
observe longer S Dor-like oscillations in this star. 



5. SUMMARY 

The aim of this paper was to investigate photo- 
metric variability of massive stars in the LMC. We 
performed a variability search am ong 599 out of 



1268 objects listed in the catalog of lBonanos et al 



(2009) which had their u nique counterparts in the 
ASAS ( Poimanski 20021 ) database, spanning 10 
years of observations. We found 13 periodic and 
104 semi- or non-periodic variables brighter than 
V = 14 mag, of which ^32% had not been listed 
as variable in either Simbad, GCVS, or ASAS cat- 
alogs. 

The majority of variable stars (73%) identified 
are RSGs and almost all RSGs exhibit some form 
of variability. The variability rate among other 
spectral types fluctuates between ~4% (for O- and 
B-type stars) and ~37% (for sgB[e] stars), but the 
absolute numbers of stars in the groups are too 
small to perform reliable statistics. 

All 13 periodic variables are EB stars, but to 
our surprise only five of them were known to be 
binary before, even though the LMC had been sur- 
veyed for bright variables about a century ago at 
Harvard. However, the sensitivity of the photo- 
graphic plates only allowed for detection of vari- 
ables with amplitudes above 0.5 mag, while all new 
EBs have amplitudes below 0.6 mag. The eight 
new EBs have O and B spectral types, which are 
consistent with their short orbital periods. We 
have already started follow-up spectroscopic ob- 
servations in order to verify their binary nature 
and investigate their physical properties. 

We present some of the more interesting 10 year 
long light curves, which for many objects are the 
most extensive light curves available in the liter- 
ature. This is particularly important since sev- 
eral of the systems are relatively rare. We also 
provide light curve stamps for all semi- and non- 
periodic variables for easy reference. In addition, 
we present tables with the stars' parameters. All 
light curves and full tables are available on request. 
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Fig. 1. — : ASAS light curves of the LBVs R 71 in V (black circles) and / (red crosses) bands, and R 127 in 
V (gray open circles), spanning 10 years. The right-hand ax is shows the absolute magnitude value, assuming 
an LMC distance modulus of 18.41 mag (|Macri et al.ll2006l ). 
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Fig. 2. — : Spatial distribution of ASAS stars in the LMC region in equatorial coordinates. Gray dots are 
the —19,000 objects brighter than V = 14 mag (out of —28,000 total) found by ASAS in this region of the 
sky. Filled circles represent the 599 final matches between the ASAS and B09 catalogs and open circles are 
the subset that was identified as variable (117 stars). Blue pentagons mark new EB candidates (see Section 
3 for details) . To show the location of the LMC bar, we also plot the stellar density contours based on the 
OGLE-III photometric maps of the LMC (Udal ski et al. 2008). We only include I < 15 mag objects, binned 
into 0.2 deg regions in both R.A. and Decl. The contour levels are 50 (for the outer) and 90 (for the inner) 
stars per bin. The outer contour includes the 30 Doradus star-forming region. Note that the OGLE fields 
do not cover the Shaple y Constellation II I star-forming region. The 1 kpc marker is based on 18.41 mag 
LMC distance modulus ([Macri et al.ll2006l ). 
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Fig. 3. — : Standard deviation cry vs. T^-band mean magnitude for all 599 matched sources. Variable stars 
are marked with open circles; EB candidates (eight points) are denoted by filled circles (see Section 3.1 for 
details). 
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Fig. 4. — : Phased ASAS l/-band light curves of the five previously known EBs, sorted by right ascension, 
spanning 10 years. The header of each light curve contains two IDs, spectral type (in brackets), and orbital 
period in days. For more information on these stars, see Table 1. 
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Fig. 5. — : Phased ASAS F-band light curves of eight new EB candidates, sorted by right ascension. The 
header of each light curve contains two IDs, spectral type (in brackets), and orbital period in days. For more 
information on these stars, see Table 2. 
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Fig. 8. — : This and the following two pages present light curve stamps of all 104 semi- or non-periodic 
variable stars identified in the course of this study. 
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Fig. 8.—: (Continued) 
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Fig. 8.—: (Continued) 




Fig. 9. — : Left panel: CMD for the 599 matches between ASAS and B09 catalogs. V and / magnitudes 
are the mean values based on ASAS light curves. Black dots represent the 117 variable stars and gray open 
squares the non- variables. EBs are marked with big open circles. Right panel: V amplitude of 117 variables 
vs. V — I color with the same designations. For the non-variables, the light curve scatter is plotted. R71 is 
marked with a cross (both panels). 
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Fig. 10. — : Infrared CMD for 113 variable stars out of 599 matches between AS AS and B09 catalogs, 
for which infrared colors were available. Values of [3.6] and [4.5] were taken from the B09 catalog. The 
designations are the same as in Figure 2 of B09. EBs are marked with big open circles and R71 is marked 
with a cross. 
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Table 1 

Parameters of Known Binary Stars 



B09 ID" 


RA 2000 


Dec 2000 


HJDo - 2450000 


Per 


Vasas 


Ampv 


Sp Type 




[dog] 


[dog] 


[days] 


[days] 


[mag] 


[mag] 




BAT99-6 (Sk -67 18) 


73.8121643 


-67.1899185 


3060.318974 


2.001248 


12.04 


0.49 


03f* + O EB 


BAT99-32 (Sk -71 21) 


80.5938721 


-71.5994720 


3425.880887 


3.815192 


12.67 


0.24 


WN6(h) SB 


Sk -67 105 


81.5257950 


-67.1824417 


3400.010490 


3.301529 


12.41 


0.43 


04 f + 06V 


Sk -67 117 


81.8641281 


-67.1984711 


3412.698954 


4.828888 


13.03 


0.66 


09 III-II + 08 V: 


MACH081.8763.8 


83.6720810 


-69.5275040 


3090.591666 


1.404787 


13.62 


0.59 


03 If* + 06:V 



"Star designations following B09: Brcvsachc r et al.l (BAT99; [19991 ) . ISanduleakl (Sk: ll970f) 



Table 2 

Parameters of Eight New Eclipsing Binaries 



B09 ID a 


RA 2000 


Doc 2000 


HJDo - 2450000 


Per 


Vasas 


Ampv 


Sp Type 




[deg] 


[deg] 


[days] 


[days] 


[mag] 


[mag] 




Sk -69 9 


72.4642105 


-69.2011948 


3445.552851 


1.736425 


12.46 


0.22 


06.5 III 


Sk -70 18 


74.0687943 


-69.9936905 


3445.140846 


7.625148 


13.04 


0.42 


B0 n 


Bl 98 


78.1266632 


-67.4539948 


3245.926372 


3.765188 


14.03 


0.62 


Bl: V 


Sk -65 47 


80.2277908 


-65.4550781 


3585.283137 


3.701956 


12.41 


0.23 


04 If 


Sk -71 29 


81.2763748 


-71.4644165 


3462.286550 


1.832248 


14.04 


0.58 


B extr 


Sk -67 102 


81.5037918 


-67.0528107 


3450.143306 


3.659535 


13.53 


0.45 


B2 III SB: 


Sk -69 194 


83.6503296 


-69.7601395 


3451.041843 


12.22083 


11.66 


0.67 


B0 I + WN 


Sk -67 270 


86.7945404 


-67.9819183 


3474.245086 


3.067998 


13.74 


0.66 


B0.5 V 



"Star designations following B09: IBrunet et~aTI (Bl: 119751 ). ISanduleakl (Sk: ll970D 



Table 3 

Parameters of 117 Variable Stars Identified in this Study 



B09 ID a 


RA 2000 


Dec 2000 


Vasas 


Ampv 


Iasas 


Amp/ 


Sp Type 


E(V-I) b 


Known? c 


Ncighb 


ors d 




[dog] 


[dog] 


[mag] 


[mag] 


[mag] 


[mag] 




[mag] 








P92 3160 


74.2139130 


-66.3970566 


12.03 


0.36 


10.35 


0.19 


early K I 






1 bright, 


~25" 


Bl 53 


74.8157120 


-69.7153549 


14.04 


0.61 






B0 III 






1 bright, 


~20" 


Sk -69 59 


75.8029175 


-69.0269470 


11.64 


0.20 






O(f) +B1 I? 


0.000 




1 bright, 


~15" 


M02 062353 


76.0413361 


-70.2049942 


12.84 


0.40 


10.86 


0.15 


Ml I 


0.120 








M02 064706 


76.2255402 


-70.5552521 


12.83 


0.40 


10.85 


0.19 


M0-1 I 


0.134 








M02 065558 


76.2917938 


-70.6675568 


12.57 


0.35 


10.62 


0.16 


M0 I 


0.000 








Sk -68 58 


77.4382935 


-68.7694168 


11.92 


0.30 






B3 lab 


0.224 




1 bright, 


~25" 


Sk -69 75 


78.3782501 


-69.5399170 


10.56 


0.51 


10.12 


0.59 


B8 I 


0.163 


V1577 






M02 119219 


80.0987091 


-69.5575867 


12.76 


1.26 


10.68 


0.91 


M3 I 


0.126 


V2148 






Sk -65 47 


80.2277908 


-65.4550781 


12.41 


0.23 


12.49 


0.27 


04 If 










BAT99-32 


80.5938721 


-71.5994720 


12.67 


0.24 


12.56 


0.37 


WN6(h) SB 


0.141 


V2296 






c GOVS LMC designation. 


d Approximate 


brightness 


and distance of close ncighb 


ors, based on DSS 


images. 


Sec Section 


3.3 for 





details. 



Note. — A full version of this table is available in the electronic edition of the Astrono7nical Journal. 

a Sta r des ignat i ons follo win g BQ9 : [Breysacher et"ld1 (BAT99; 1 19991) [Sanduleakl ( Sk; 1 19701) . [Wester lundl ("W; [TMTl ) . IBrunet et~aT| (Bl: 119751) , 
ILuckd (LH; fl97^) , |Henizel (S: fl956l ) . [Parker et all fP92; |200l|) . IMassevl fM02; 12003) 

b Mcan extinction values taken from Pc.icha & Stanck (2009). Zero values (0.000) correspond to negative values for the calculated extinction 
A/, which is the case for 30 out of 101 stars for which the extinction value was available. 
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